















章数学预备知识 



































证毕. 



















































要有效地求多项式的值，我们可以用嵌 if 乘法把^分组： 

































































































































































































































































































































































































非线性方程的解 





























































































































































































































































































这个方程组的解就是 







-对非线性方 


第3聿 



H7 以(《4, v ) = (.3, 一 4) 为初值，使用贝尔斯托法求»例中的多项式的实二次因式. 
* 以双精度形式编写基于上面伪代码的计算机程序，并且在第7次迭代，得到 
















我们有 















































































































4 章解线性方程组 


































































为了解这个方程组，傾定对一切 i, a ^ O , 然后，从第一个方程得到: c,. 把巳知的这个 or ，值 
代人第二个方程，得出第二个方程的解用同样的方法依这个次序一次得到一个: t., 最后 
可得到 x,, x 2 ，… ，： 此时求解的一个有效的算法称为肉餺回代 ■ 


习惯上，在#< 0 的情况下，形如的任何和式被认为是 0. 

同样的思想可用来求解具有 ik 角鏑 掬的方 程组. 这样的矩阵方程组具有形式 




必须再次 餒定： 对 a .^ 0 . 求解： r 的一个有*的算法如下，并且称之为 ft 后回代. 


Mtpat (x,> 

还存在另一种 


种利用上述思想容易求解的简单的方程组类型- 
_的方程所得到的方程组.为了说明这点，考察方程组 
fan a it 0 I 


只需对这些方程重新排序，我们就可得到一个 T 


-即，置換一个 H 角方程组中 












在 ）</, 且对一切 f, 的情况下，这个算法起作用. 

上述 W 个算法都没有给出直接转换成大多数程序设计语 豸所霱 的足移细节. 
三角方程组的向后回代算法可能就需要更为详尽的指令： 








































































































































































































































































































































































-般公式是 



























































解线性方程钮 




















































































































































数值线性代数精选 






特征方程的根是 


•(A-4>a*+A + 2) =0 






















因为等式 ⑸ 中分子收敛于 0 而分母收敛于 ( C -l)* 且不等于 0, 所以显然有 r>/A«=0. ■ 
因为公式中减法相消将最终损害结果，所以重要的是在明显地出现稳^值后，立刻终止艾 





































1 舒尔分解 
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在这两个矩阵中，未显示的元素都为 0. 通过首先作>1的一个奇异值分解 
A = PDQ 

再取 

a + = Q'ir p* 

来定义一个一般矩阵 A 的广 义逆. 在后面我们将看到虽然一个矩阵的奇异值分解不唯一，但 
是其广义逆却是 唯一碥 定的. 

例3例1中矩阵 A 的广义逆是什么？ 

解因为此例中 P 和<?是单位阵，由此直 接耐得 



室成下列矩阵形式的方程组 


I小欧几里得范败的元素. 


I小二# 




S 中具有最小欧几里得范败的 


另一种定义叙述 如下： 利用欧几里得范败，设 

/ > = in/< || Ar-*|| 2 « x6 C-} 

则方程的 最小解是集合 K = U» II Aa ■—6 沁中最小欧几里得范败的元素.显然易 
见这个定义包含了前面所描述的四种情况.例如，若 #> =0,我们有情况1和2,而情况3和4 
別对应于彳>0. 

定理 2( 广义逆最小解 定理） 




























































































































第 6 幸 


(回忆宪罗内克 S 函败的定 义是： 当 = « 时心 =1，当*乒£时夂 =0.) 我们不难得到一组具有 


…， x. 的集合， 这些多项式被称为基函 ft . 而用它们， （9) 式给出了拉格朗 

中数据的基函数以及拉格朗日 a _ 值多项式是什么？ 

15, -7, -6, 0. 所以基函敫是 

ri= (^7Mf+«) > r-so x(i+<i)( - | - 7) 


拉格朗 H 型插值多项式是 

pj < x ) =&( x )- 23/, < x ) - 54/, ( x )- 954/, ( x ) 

这个表达式看上去与 (7) 式中的 p , 不同 • 尽管它们的表示形式不同，但作为一个*教来说它们 
是恒等的. ■ 

例2求出下列两点败表中的插值公式. 



























































































函教遥近 











: 1 丄 | 4此 






















































































































































































































A 數逼近 


a 公式成立.在此假设之下，我们来《实下一个情况.归纳鋇设的是 do) 式，利用前面所说明 
的记号把它重写为紧凑形式： 

£ b * = iatBr (11> 

根据基本递归关系 .(3) 式，我们有 

差 sr = £[vr>B?+a-va' )^,3 

=vr £b ； + (1 - V!t!) £bI, -attlBl, 02) 

其中两次用到了 （7) 式. 接下来，我们归纳飯设 （11) 式替換出现在 （12) 式右蝴的导数.其结 
果是 

»vr ! (fc?Br — bw ) +bj 
+ (1 - VJf; HhaUB^l - haUzKD-aSlBl, 

通过简单的重新组合，上式可以写成下列形式 

=«r'B ? + taivrBr -«&'%, _ 

- fc^a-VS,' fe WBW < 13 > 

+ fe^,d-vw)iKr 1 , 

我们将对这个等式做下列替换，那么所有的证明都可以利用 (8) 式和 (9) 式来完成， 

hivr-'Br =fcr'v}Bf-' 

- fe^d-VS,' )BK =— tejt/d- V^)BS,' 

- k „ U , VT ' Bi ；}+ kaU , (1 _ VJf/ )Bt； 

= kal,(l -WBW 

现在 （13) 式变换后的形式为， 

£fir' =«r'Bj+*[«^*vjBr> +«j«(i-vj„)bw] 

-*[„{；/ VJhBJt / +^(1-^,)^*] (14) 

利用基本的通归关系 (3) 式，化简 （14) 式带括号的項，得到 

£Br - 瑞 B* m - 

=« + DarBi - a + (15) 

因为这个等式是用 a +1 替换了 （m 式中的 a , 所以归纳证明步驟完成. ■ 

引理 6(B 样条光滑性51 理） 对于 B 样条 B? 属于连续函教类 C*-i(R>. 































































































































































































































































































•t.)=«+00 而且茗 （i,>=/(JE,). 如果当 J： 离开: C, 时 
S 对 g(X,> 几乎没有彩畹. 

的一 种选株 ft 




如果= llx - x . ll 2 ，那么去掉奇点以后，叫有基性质 :它在点：， 取值为！而在其 
他结点处取值为 0. 

6.10.10 多霣二次插值 

另一个多 元插值过程是 R . L . Hardy 提出的.它的基函数是所谓的多重二次 * 数 

«.(/•) = { II /•-/»< II*+c*>« rt (1<,<») 

其中的范数是欧几里得范数， r 是参数， Hardy 建议 c 等于结点间平均炬离的 0 . 8 倍在利用 
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6.12. 3内积，伪内积，伪范数 

复希尔伯特空间1^[—中的内积定义为 

由(其中 A=0, 士 1, 士2,…)定义的函数 E, 构成一个复希尔伯特空间中函数的标 
准正交系.这意味着当时 <&• 艮>=0并且 <E*, £,> = 1. 这可以从下面的计算中看出. 
当㈣时 

<E k ,E,>= ij* E t (x') E a (x)dz = 去 J e a *e'* ,I dr 



显然地<£,， E 4 > = 1. 而且函数 &构成 一个标准正交序列. 

使用下列内积符号将会很方便 

</.«>« - jmjm (12) 

因为从条件</, /> N = 0 不能推得/=0,只能得知在毎个结点 2W/N 上两数 /(x) 的值是0, 

所以这个函数不是一个真正的 内积. 但是它符合(复）内积的其他性质.它们是， 

复内积性质 

1. </. f>N>0. 

2. </. />!». 

3- ( af + pg , h ) s - a < f , h ) N + p < g , h > M . _ 

与伪内积一起，还有下面定义的伪范数： 

II /In =■✓</./>« 

我们有II / II N = 0 当且仅当 f (2 nj / N )~0 . 其中 
对于插值而言，卜述定理起着决定性的作用. 

定理 2( 伪内积定《>对任意我们有 

, P P , /I 若/V 整 

| o 其他 (13) 

证明问埋中的表达式 町写为 

如果/V整除*_»«，那么是整数并且 e w-"N = l. 因此，每个被加数是1，它们的 
平均数也是 1. 另一方面，如果不能被 iV 整除，那么这时我们可以应用 
几何级败求和的标准 公式： 

其结果是 

































































































































































































































































































狭出收轚性. 


自适应求积方法是指自动地利用被积函数的性质来计算定积分，理想情况下， : 


用户只需提供被积函败/、区间 [ a , 6] 以及所要求的精度 e . 然后程序将区间划分为各种长度 
的小段使得在子区间上数值积分产生満足精度要求的结果. 

这里讨论一个典型的自适应求积方法，主要用到子区间上的辛普森法則连同所涉及的误差 
估计. 7 . 2节中（6>式给出辛普* 法南： 


£/(x)4t= S(a,W- ^[(6 —a 

S(a ,»= ^[/(a) + 4/(^) + /(*)] 












































































































































































根据微分方程，我们有 

xiO= / 

x'U') = /, + /* X ’ =/,+/!/ 

x *0)= /. +/«/ + (/«+/,/)/, +/(/-+/-/> 

这里下标表示偏导数，反复地利用微分法的链式法 M. (2>式中的前=.项现在可写成下列形式 
xO + fc) - x + A/ + (/, + //, ) + 0 (A*) 

-j： + i-A/ + |A[/+A/. + A//x] + O(A') ⑶ 


其中 a ： 表示 x (»)，/ 表示/0, *>, 等等.我们借助于两个变童的泰勒级数（见 
几项可消去偏 导败， 

fu+k,x+hp = f+hf.+hff.+oa*、 

(3〉式可改写成 


x(< + A> =x + jA/ + + A,x + A/) + CKh” 

因此 I 步进求解的公式是 

x(/ + A) ■= x<*) + + +A,x + ) 

或等价地， 

x(< + A) = x(t) + j(F,+F,) 


l 节）中的前 


(4) 


|F, ■ A/(i,x) 國 

If* = a/(<+a,x+f,> 

这个公式可反复地使用，每次前进一步 求解. 它被称为二阶龙格-库*方法，也称为 H«m 
方法. 

一般说来，龙格-库塔公式具有形式 

«(* + A ) - x + to , A / + WthfU + Ji,x+phf)+Oa*) (5) 


















C 是一个与 ft 
C 不变.设 T； 



由下列形式的公式 算出： 




















































































现在考虑初值问埋 








































































































































常歡分方程軚值解 




c 复变* 2 导出的 • 为了现察〆的级数的收敛性，我们在 c •上取 任意范数并使用 




当 || A || 时.这个最后的表达式是通常的指数级数的尾部.因此，当 woo 时，〆的级败 
的尾部收敛于 0. (众所周知.这个理由采用了给定范数的”X"矩阵空间的完备性 .） 

若* 是一个实变量，則由我们的定义得出 




































































































9 章偏微分方程数值解 








































































































达到抜小： 













































































































































































































































+ - 2 v h +i 4 . y .,) = A*g f h aiaik * 

由引理 1 知，具有元索 siniV 的矩阵非竒异.所以，可以从(7>^^出 

^ (- 4sin J ^ )+ , - 2^ + i*,,- , = A* /„ (8> 

初看起来 (8) 式好像是另一个有 „* 未知量和 ” 2 个方程的方程组，它仅仅稍稍不间于原来的方 
程组 (3). 但仔细的检査 显示： 在 (8) 式中， A 可取固定的，并且因为所得的《个方程的方程组 
是三对角的，所以它容易直接 求解. 因此，对固定的 h (8) 式中的未知量构成 R- 中的一个 


上面所用的方法把原来的” 2 个方程的方程组拆开成„个方程组，而每个方程组有„个方程 . n 
个方程的三对角方程组可用 0OO 次运算求解 （事 实上，需要少于10«次运算).因此，我们可 











































































































































































































































































































































附录 A 数学软件一览 
































































































































































这是求解已知初值的常微分方程的子程序 集成. 这些程序对用线方法求解不定常偏微分方 
程时形成的常微分方程的求解特别有用.它们是由 Lawrence Livermore 国家实验室 （LLNL) 所 
开发的，这个程序包的主页是 

http,//www. Uni. gov/CASC/odepack 

A. 5.6 PDE2D 

这是在一般的二维区域和三维盒体内 求解线 性和非线性不定常、定态及特征值偏微分方程 

组的一个有限元代码.它能处理诸如掸性、扩散、热传导、势能和流体力学等领域中的 问®. 

它具有一个交互接口和大量的图像输出能力.它在 

http；//members, aol. com/pde2d/ 

上可以得到. 


这是求解偏微分方程及相关问理的科学计算的一个便携式可扩充的工具箱.它是一套创建 
大范围应用的败据结构和程序，这些应用由并行的线性和非线性方程求解程序和无约束最小化 
模块所 组成. Argonne 国家实验室 (ANL) 支持这个项目并且是主要的发布站点， 
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